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AN IMPROVED METHOD FOR CALCULATING SUPERSONIC 
PRESSURE FIELDS ABOUT BODIES OF REVOLUTION 

By Robert J. Mack 
Langley Research Center 

SUMMARY 

An improved near- field method for determining supersonic- flow- field proper- 
ties about a body of revolution is presented and discussed. Comparisons between the 
improved method, Whitham's theory, and wind-tunnel results are shown for four bodies 
of revolution — three closed-nose bodies and one ducted body. At Mach numbers of 2.96, 
3.83, and 4.63 and ratios of radial distance to body length of 1.0, 2.0, and 5.0, results show 
that the improved method does reasonably well in predicting flow-field pressure signa- 
tures and represents a definite improvement over existing near-field theory. It is also 
shown that the simple ’’area balancing" shock prediction technique works reasonably well 
for bodies with pointed nose sections provided the improved method is used to locate 
the body flow-field disturbances more exactly. 

INTRODUCTION 

Supersonic- flow- field prediction techniques based on the modified linearized theory 
of G. B. Whitham (ref. 1) have proved to be successful in the low supersonic to Mach 2.0 
speed range. Reference 2 demonstrates this capability for a variety of bodies of revolu- 
tion and reference 3 shows it extended to aircraft and wing-body combinations. 

Recently, wind-tunnel tests were conducted to explore the applicability of existing 
theory in the high supersonic, low hypersonic speed range (ref. 4). Near-field theory 
pressure signatures compared well with the measured pressure signatures in the wind- 
tunnel tests at a Mach number of 2.96. However, disparities between the theoretical sig- 
natures and the wind-tunnel measured signatures appeared in the tests at Mach 3.83 and 
became more pronounced in the tests at Mach 4.63, especially for the low- fineness- ratio 
bodies and the blunt-nosed bodies. The lack of agreement occurred in nose-shock loca- 
tion, signature shape (particularly pressure gradients), and integrated area under the com- 
pression part of the pressure signature. 

In this report, corrections to the near- field Whitham solution are presented. Com- 
parisons of wind-tunnel data, near-field theory, and the improved near-field method are 
also shown. Although these linearized theory solutions are strictly applicable to pointed 



bodies (A’(0) = 0) or slender ducted bodies (R(0) ^ 0 and R'(0) small enough to generate 
linearized theory disturbances), a method is outlined for use on blunt-nosed bodies 
jA’(O) - <» or R’(0) S ~j. 

SYMBOLS 


A body cross-sectional area 

Cp pressure coefficient 

f(t) source distribution function 

F(y) Whitham's F function (ref. 1) 

F(y) impulse function defined by equation (8) 

h area influence function (called h(x) in ref. 1) 

I: _ (r + 1)M 4 



V-ML 

XV " i 

l reference body length 

^ ducted-nacelle length 

L(y) initial source strength function defined by equation (9) 

M Mach number 

p free-stream static pressure 

Ap incremental pressure due to flow field of model 

Ap max maximum value of Ap 
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r 

R 

Ri 

R f 


Rmax 


R, 



y 


y Q 

y l’ y 2 


radial distance from body axis 
body radius 

body radius at initial Mach plane intersection (see fig. 3(a)) 

body radius at final Mach plane intersection (see fig. 3(a)) 

maximum body radius 

radius where Mach cone is tangent to body 

dummy variable of integration in x- direction 

distance in longitudinal or windward direction, measured from model nose 
distance x to where Mach cone is tangent to body (see fig. 3) 
distance x to tangent Mach cone apex (see fig. 3) 

(see fig. 3) 

value of characteristic defined in x-direction on body by y = x - /3R(x), and 
in flow field by equations (1) and (7) 

value of y for which ^F(y)dy is maximum 

values of y which uniquely satisfy the trailing- shock equations 


0 = /m 2 - 1 


y ratio of specific heats (1.4 for air) 

M Mach angle 

A prime is used to indicate a first derivative. The symbol A denotes incremental 
quantities. 


3 



THEORY 


From reference 1, the equation of the characteristic (i.e., the loci of points influ- 
enced by a given source distribution) can be written as 


= x . fr + k±JJMiJ y 




Vy - t + 2/3r - jy - t + 2/3R(y) 




f’(t) dt 


2M 2 ^ 


y ln Jl Vy - t + 2j3r - h - t 


^y - t + 2|3r + /y - t 


l/y ~ t + 2/3R(y) + ^y - t 


[l/y - t + 2£R(y) - 

Using the assumption that pr/y is large, equation (1) reduces to 
y = x - /3r + kF(y) ft 

where 


. f T (t) dt 


F(y) = ^ 

J 0 ][■ 


(t) dt 


i/F~t 


In this paper, the Stieltjes integral form 


F(y) = r y r_2_1 1/2 
w Jo \^R(t )) 27T 


(l) 


( 2 ) 


(3) 


(4) 


developed in reference 1 for bodies with discontinuous surface slopes is used. 

The nose shock angle is determined from the assumption that, to first-order cer- 
tainty, the shock bisects the free-stream characteristics and the body surface character- 
istics. With the use of equation (2), reference 1 gives 

kF 2 (y) |/r = 2 f Y F(y) dy (5) 

J 0 

as the condition satisfied along the nose shock. A similar treatment, also developed in 
reference 1, results in 
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as being the conditions which determine the positions of the trailing shocks. 

Complete details of this "area balancing" technique are found in references 1 
and 5. Notice that in equation (5) the difference between the value of the characteristic 
and the distance that the shock stands ahead of the Mach line from the nose (that is, 
y - (x - /3r)) multiplied by the function F(y) sets the value of the right-hand side of the 
equation which in this case is simply twice the integral of the F(y) function. This par- 
ticular feature also shows up in the improved equation for the location of the nose shock. 
Notice also that equations (2), (5), and (6) express the characteristic equation and the 
shock conditions as being functions primarily of distance and secondarily of body shape 
via the function F(y). 

An approximate near-field equation of the characteristic is obtained by applying the 
method of integration by parts to equation (1) in order to get an expression containing the 
function F(y) and then by using some empirical simplifications to reduce that complex 
relationship even further. The equation that results from these operations is 

y ~ x - /3r + kF(y) [r 1 / 2 - R 1//2 (y)] + EF(y) [r- 1 / 2 (y) - r - !/ 2 ] + L(y) (7) 


where 


F(y) = = 2f(0) /y + f y F(y) dy 

J o J o 


( 8 ) 


and 


L(y) = 2M 2 f(0) In 


fr i 

n -il 

/y + 2/3 r + tfy 

l/y + 2/3R(y) - /y 

\p + 2|3r - \Jy 

Vy + 2/3R(y) + /yjj 


(9) 


These expressions show that both the body radius R(y) and the integral of the func- 
tion F(y) are significant factors in the near-field; R(y) is the radius where a rear- 
ward running Mach line from y intersects the body surface and is found analytically or 
by numerical iteration from 


y = x - /3R(x) 
R(y) = R(x) 


( 10 ) 


Note also that the Whitham boundary condition y = x-0r at r = R from refer- 
ence 1 is satisfied by equation (7). 

The integral of the function F(y) appears in equation (30) of reference 6 where it is 
seen to influence flow-field predictions even at large radial distances. It also is found in 
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a somewhat modified form in reference 7, which gives a method for extrapolating near- 
field pressure signatures to the mid- and far-field. 

Numerical tests were performed to determine how well equation (7) approximated 
the behavior of equation (1) in the near-field, with a 5° cone taken as the test case. The 
results showed that through a Mach number range of 1.25 to 5.0 and a range of r/x from 
0.09 to 0.17, equation (7) was sufficiently accurate to warrant the work necessary to 
develop a set of shock equations. 

Using equation (7) and the shock bisection assumption gives 
kF 2 (y)[r 1 / 2 - R 1 / 2 (y)] + KF(y) F(y)[ir 1 / 2 (y) - r -1 / 2 } + F(y) L(y) 


= 2 F(y) dy + | f Y F 2 (y) R'(y) R _1 / 2 (y) dy 
J 0 l Jq 

+ kJ o y [R- 1 / 2 (y) - r-l/ 2 ][F’(y) F(y) - F(y) F'(y)] 

+ 1 f Y F(y) F(y) R’(y) R" 3/2 (y) dy + f Y [F’(y) L(y) - F(y) L’(y)] dy (11) 

2 J 0 j 0 

as the condition satisfied on the nose shock. Although complicated, equation (11) can be 
programed on a digital computer so that nose-shock solutions can be easily found. 

A similar treatment with the intersecting characteristics originating on the body or 
wake surface and with some simplifying assumptions gives 



|(y 2 - yi)[ F (y2) + F (yi )] + 1 F (yi) F (y2)[ Rl/ 2 (y2) - Rl/2 
+ |^ 2 ) R ' 1 / 2 (y 2 ) + F (yi) R " 1 / 2 (yi)][ F (y2) - F (yi)_ 


i^i) F H - F (yi)] 
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and 


kF 


( y i) 


= kF/ 


*1/2 - R 1 / 2 ^ + kFjy^ R" 1 / 2 ^ - r-1/2] + L( Yl ) - yj 
(y 2 ) rl//2 - Rl / 2 (y2) + kF(y 2 ) R_1 / 2 (y 2 j - r" 1 / 2 ] + L^y 2 j - y 2 (13) 


as the governing conditions on the trailing shocks. An examination of equation (12) shows 
that "area balancing" is the core of the solution; area balancing for the improved method 
is represented by equation (11), with only the first term on the right-hand side being kept 
and used, and by equation (12) reduced to 

Q F « dy 4 ( y 2 - y l) If fa) * F ( y l)] 

The other terms in equations (11) and (12) are combinations of body radius, surface 
slope, F(y), and free-stream Mach number which serve to shift the shock location for- 
ward or rearward from the location established by the simple area balancing solution. 

A comparison of the area balancing method and the shock equation method is made at the 
end of the next section. 


RESULTS AND COMPARISONS 

Four bodies of revolution are used to compare wind-tunnel data, near- field theory, 
and the improved near-field method: a cone-cylinder model, a model with a linear-area 
forebody and a cylindrical afterbody, a model with a "lower bound" (ref. 8) forebody and a 
cylindrical afterbody, and a truncated cone- cylinder ducted-nacelle model. Although the 
wind-tunnel models are of finite size, semi-infinite mathematical representations are used 
in the computer program. The pressure disturbances from the model sting are excluded 
from the signatures, so the wind-tunnel models are, in effect, also semi-infinite. 

The flow- field pressures are calculated from 

Ap = yM 2 F(y) ( 14 ) 

P ^2/3r 

(which was obtained from ref. 1) so that a direct comparison between near-field theory and 
the improved near-field method can be made. Since most of the data are obtained from a 
sonic-boom report (ref. 4), the measured signatures and the near -field theory signatures 
are shifted so that they are alined with the improved-method signatures in the same 
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manner and for the same reason as in reference 4. Also, similar dimensionless pressure 
and length parameters are employed to avoid confusion and aid direct comparison. 

The following order is used to present and discuss the models, test data, and theo- 
retical predictions: (1) cone- cylinder model, (2) linear-area and "lower bound" models, 
(3) ducted-nacelle model, and (4) comparison of the area balancing and the shock equation 
solutions of nose-shock location. 


Cone- Cylinder Model 

Flow-field pressures are calculated for the lowest-fineness-ratio cone-cylinder 
body in reference 4 (see model 1 in fig. 1) and are presented in dimensionless parameter 
form in figure 2. This body, for which the ratio of the semivertex angle to the Mach angle 
is the largest in the cone- cylinder series at all Mach numbers, constitutes the most 
severe test of the improved near-field method. Mach numbers of 2.96, 3.83, and 4.63 and 
values of r /l of 2.0 and 5.0 are used to compare the experimental data and theoretical 
predictions (see fig. 2). No boundary- layer displacement thickness is added to the body 
radius because a complete and accurate knowledge of the flow over this as well as the 
other bodies of revolution used in this report is not available. 

At Mach numbers of 2.96 and 3.83 with j = 5.0, the predictions from the improved 
near-field method compare well with wind-tunnel data (see figs. 2(a) and 2(b)). Peak pre- 
dicted overpressures are higher than measured overpressures possibly because of the 
rounding effects of wind-tunnel turbulence and model vibration, but nose-shock location 
and signature shape show the differences caused by the extra terms in the equation of the 
characteristic (eq. (7)). 

Agreement between experiment and prediction is less complete at M = 4.63 with 
r/l values of 5.0 and 2.0 (see figs. 2(c) and 2(d)). At this Mach number the ratio of cone 
semivertex angle to Mach angle is 0.516; this value is above what is considered the limit 
of linearized theory. Nevertheless, the improved near-field method does provide the 
better solution and picture of what is happening in the flow field. Particularly interesting 
is the fact that no distinct trailing shock is measured or predicted by the improved method 
at the test Mach numbers and at the distances selected for measurement. This feature is 
probably due to the slow coalescence of the body characteristics in the region of the shoul- 
der plus the tendency for wind-tunnel turbulence and model vibration to smear out the 
recording of any small shocks that might have formed. Signatures with nose shocks but 
without trailing shocks are also presented for the other models in this report. 

Linear-Area and "Lower Bound" Models 

In addition to the cone- cylinder model, figure 1 also shows the next two models (from j 

ref. 4) which were used for flow-field comparisons. Like model 1, models 2 and 3 have ; 

j 
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the lowest fineness ratio in their respective series in reference 4. Since they are blunt 
^R’(O) ^ , a somewhat different treatment must be used to get flow- field information. 


When a blunt body is "Mach sliced" in order to determine its disturbance-oriented 
area distribution and therefore its source strength function, it is found that the sources 
begin at an initial value of x Q = x^ - /3R^, as shown in figure 3. As the Mach planes move 
rearward between x = x 0 and x = 0, the areas intercepted are entirely on one side of the 
body axis. Body radii Rf and Rj are the maximum and minimum values on the merid- 
ian section shown in figure 3(a). However, the derivation of equation (4) assumes that 
only one intersection exists, so the nose is suitably modified to make the theory applicable 
over the entire length of the body. 


The method that gives good results, that is consistent over the Mach number range 
described, and that is plausible in concept is shown in figure 3(b). A tangent cone with its 
apex at x^. = x Q - x^ = -0R^ is used to start the new representative body of revolution. 
Figure 3(c) shows the components of the F(y) calculation and how those of the represen- 
tative body (blunt body with tangent cone) compare with corresponding values of the blunt 
body. 

A point in the flow field cannot "feel" the influence along Mach or characteristic 
lines of a disturbance from the nose of the blunt body because of the shielding that is pro- 
vided by the nose itself. Any initial disturbance experienced appears to originate from 
sources ahead of the nose. Therefore, putting sources (body volume) in a tangent cone 
does not change the flow- field representation if the new function F(y) has approximately 
the same shape and, more importantly, the same integrated area as the blunt body, as is 
shown in figure 3(c). 

This tangent-cone method insures that R’(y) is finite everywhere, that R(y) has 
a single value everywhere, and that the distance to the tangent cone apex x Q is fixed by 
body geometry and Mach number considerations only. The results from using this method 
on models 2 and 3 are shown in figures 4 and 5. 


It can be seen that the wind-tunnel data and the theoretical predictions compare less 
favorably as bluntness and Mach number increase, but that for a body of moderate blunt- 
ness (model 2) the improved near- field method does reasonably well over the Mach num- 
ber range covered by the tests. 

A convergence test of the tangent- cone method (based on varying A l/l) is shown in 
figure 6. Model 3 was used, with M = 4.63 and -f = 5.0 as the test conditions. The 

A 7 ^ 

reference line of 1.0 is set for -j- = 0.02 because the number of area stations used to 
describe the forebody is large (51) but not excessive. 

The lines drawn between the circled points in figure 6 are trend indicators rather 
than loci of values and are used to show that for this blunt model, Mach number, and radial 
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distance, a very large number of area stations are needed to reach convergence. Most of 
the scatter in the maximum-pressure- ratio data points is due to the effect of small differ- 
ences in the individual F(y) functions. The steep pressure gradient following the nose 
shock coupled with the numerical technique used in solving the nose-shock equation mag- 
nifies these differences and causes the observed behavior. However, the integral of the 
positive lobe of the F(y) function (the shaded area shown in the inset sketch) is much 
more stable. When the pressure signatures for which ^ fk 0.01 are overlaid, scarcely 

A l t 

any difference can be seen. If the one data point at -y = 0.01333 is excluded, the use- 
ful range could be extended to 4^ = 0.02. Thus, while a value of 4^ = 0.004 is needed 

L A7 l 

to produce the most accurate solution, — = 0.02 gives a satisfactory and usable one. 

a ; t 

For this reason, a value of -y of 0.02 is used to obtain pressure signatures for both 
blunt bodies in this report. 


Ducted-Nacelle Model 

Figure 7 shows the ducted-nacelle model used in reference 9. The equation of the 
characteristic used in that reference and in that application was 

y = x - /3r + kF(y)jr-*"/^ - R-*-/^(0)] (15) 

The area balancing technique of references 1,5, and 9 was used to locate the shocks in the 
flow field surrounding the nacelle. 

It is interesting that both the method used in reference 9 and the improved near- 
field method do equally well in predicting nose-shock location and strength, as is seen in 
figure 8. However, the pressure gradient in the flow field caused by the shoulder on the 
nacelle is better predicted by the improved method, even at M = 2.96. 

Comparison of Shock- Location Methods 

Figure 9 shows a comparison of the nose-shock locations and strengths calculated 
by the area balancing method and by the shock equation method. The same form of the 
characteristic (eq. (7)) is used to predict the nonshock part of the pressure signature. 

With the exceptions of the very blunt model 3 and the ducted-nacelle model (see sig- 
natures in figs. 9(c) and 9(d)), there seems to be little justification in favoring one method 
over the other. However, these exceptions point out the factors which alter the area 
balancing prediction. 

The nose shock on the ducted-nacelle model starts out at the lip strongly two dimen- 
sional in angle and strength, whereas on model 3 it is theoretically a strong detached bow 
shock. In both these cases, the effects of radius, surface slope, and fineness ratio are 


10 



very pronounced. Therefore, it would seem that the area balancing method is applicable 
only to bodies of revolution with conical nose sections. 

CONCLUDING REMARKS 

An improved near-field method for calculating the pressure disturbances in the 
supersonic flow field about a body of revolution has been presented and discussed. Com- 
parisons of the predictions of the improved method with wind-tunnel data and Whitham's 
theory have been made with the use of three closed-nose bodies of revolution and one 
ducted body of revolution at Mach numbers of 2.96, 3.83, and 4.63. It has been demon- 
strated that reasonably good flow-field predictions can be made with this improved 
method; thus the increased complexity in the system of equations used to predict the pres- 
sure signatures by this method is justified. Although emphasis has been placed on sonic 
boom (because experimental data were so readily available), the improved method is 
applicable to the calculation of near-field supersonic-flow properties about pointed or 
ducted bodies of revolution. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 24, 1971. 
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Figure 2.- Pressure distributions in the flow field of model 1. (Experimental data taken from ref. 4.) 





















(c) Sketch of F(y) calculation factors. 
Figure 3.- Concluded. 
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Figure 4.- Pressure distributions in the flow field of model 2. 
(Experimental data taken from ref. 4.) 
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Figure 4.- Continued. 
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(a) Model 1. 

Figure 9.- Comparison of theoretical nose-shock location and strength. 
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(b) Model 2. 
Figure 9.- Continued. 
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(c) Model 3. 
Figure 9.- Continued. 
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